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RESEARCH MEMORANDUM 

EFFECTS OF JET EXHAUSTS ON FLIGHT -DETERMINED 
LONGITUDINAL AND LATERAL DYNAMIC STABILITY 
CHARACTERISTICS OF THE DOUGLAS D- 558 -II 
RESEARCH AIRPLANE 

By Chester H. Wolowicz and Herman A. Rediess 

SUMMARY 


A flight investigation using pulse techniques has been made to 
determine longitudinal and lateral dynamic stability characteristics of 
the D- 558 -II research airplane with particular reference to the jet 
exhaust effects of the rocket engine. 

Any cylinder firing combination tested that included the top cylinder 
resulted in a comparable loss in directional stability. These effects 
were most severe at the highest test Mach number of approximately 1.6. 

With only the two middle cylinders firing (horizontal plane), the power 
effects were small. At the higher supersonic Mach numbers the large 
adverse power effects on directional stability were insensitive to pres- 
sure ratio from 4.4 to 15*2 in the Mach number range from 1.35 to 1*56. 

The power effects cause the rudder to float into the relative wind 
during power-on yawed flight and, as a result, tend to offset the desta- 
bilizing influence of the jet exhausts evidenced during rudder -fixed, 
yawed flight. 

The results show that the longitudinal stability of the airplane is 
not influenced by the jet exhausts, probably because the horizontal tail 
is outside the field of action of the jet exhaust effects. 


INTRODUCTION 


During the early phases of the stability and control flight test 
investigations of the Douglas D- 558 -H research airplane, it was discov- 
ered that the jet exhausts of the rocket engine had an influence on the 
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rudder hinge moments at supersonic Mach numbers (ref. l) . A general 
experimental study of this influence (ref. 2) concluded that existing 
strong external shock waves at the jet exit during power-on conditions 
caused asymmetric flow fields at the rear of the vertical stabilizer. 

The results of a subsequent wind-tunnel investigation of the effects of 
a cold-jet exhaust upon the flow over the vertical stabilizer of the 
D-558-II airplane model at several supersonic Mach numbers showed a small 
influence on the lateral stability derivatives C v , C~ , and C 7 

x p n p 

(ref. 3). The behavior of the actual airplane, however, suggested the 
possibility of more pronounced power effects. As a result, particular 
emphasis was placed on the power effects in the flight investigation to 
determine the significance of this influence on the airplane’s stability. 

This paper presents the results of the aforementioned flight inves- 
tigation. Emphasis is placed on the effect of jet exhaust on the longi- 
tudinal and lateral stability characteristics of the airplane over a 
Mach number range from O.67 to 1.6l. The supersonic longitudinal and 
lateral data and the subsonic lateral data were obtained with the 
D-558-II (l44) rocket -powered airplane. Because of relatively little 
subsonic longitudinal test data available for this airplane, the data 
from the D-558-II (145) airplane - powered with a turbojet engine as 
well as a rocket engine - were used to provide a more complete coverage 
of the Mach number range. The data were obtained at altitudes of 20,000, 
30,000, 45,000, and 60,000 feet at a nominal lg load factor. 

The lateral flight -determined power-off stability derivatives are 
compared with available wind-tunnel derivatives in the appendix. In 
addition, the lateral derivatives are also compared with calculated 
derivatives . 


The results of this investigation are referred to the body system 
of orthogonal axes. 


SYMBOLS AND COEFFICIENTS 


normal acceleration, g units 


a t 


corrected transverse acceleration, g units 


a. 


indicated transverse acceleration uncorrected for 
instrument position, g units 


b 


wing span, ft 



rudder hinge -moment coefficient, H L n . S e momen ^ 


2qMr 
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dC 


rudder hinge -moment parameter. 


hr 


dp 


lift coefficient. 


Lift 

qS 


a 


lift coefficient for trim 1 g flight, 

qS 

dC T 

lift -curve slope, — — , per deg 


da 


rolling -moment coefficient. 


Rolling moment 


qSb 


dC, 


damping -in -roll derivative, , per radian 

dE^ 

2V 


C *r 


"m 


rate of change of rolling -moment coefficient with yawing 

dC 7 

angular velocity fact or , — per radian 

2V 

dC, 

effective dihedral derivative, , per radian 

dp 

pitching -moment coefficient, f 2 

qSc 


dC_ 

, per radian 

" dS 

2V 


“q AqC 


dC. 


longitudinal stability derivative, — — , per deg 

da 


■'m 


a 


dC m 

C m, = — — , P er radian 
u ^dc 


2V 


yawing -moment coefficient, 


Yawing moment 
qSb 




rate of change of yawing-moment coefficient with rolling 

dC n 

angular velocity factor, — -, per radian 

2V 
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C 


Y 



C l/2 


I 


Y 


I XZ 


I 


z 


V Iz o 


rate of change of yawing-moment coefficient with yawing 

dCv, 

angular velocity factor, — -, per radian 

3Tb 

2V 

dCr, 

directional stability derivatives, — — , per radian 

ap 


rate of change of yawing -moment coefficient with rate of 

dc 

change of angle -of -sideslip factor, — per radian 

2V 


lateral-force coefficient. Lateral force 

5S 

dCy 

lateral-force derivative, — ±-, per radian 

dp 


number cycles to damp to one-half amplitude 


mean aerodynamic chord, ft 
rudder pedal force, lb 
acceleration of gravity, ft /sec 2 
pressure altitude, ft 

moment of inertia of airplane about body X-axis, 

P P P 

IXq cos^e + Izq sin^e, slug-ft^ 

moment of inertia of airplane about body Y-axis, slug-ft 2 

product of inertia referred to body X- and Z-axes, 

- X X 0 ) sin 2e, slug-ft 2 

moment of inertia of airplane about body Z-axis, 

2 2 2 

Iz cos e + sin e, slug-ft 

moments of inertia of airplane about principal longitudinal 
and vertical axes, respectively, slug-ft 2 


M 


Mach number 
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M r 

m 

P 

P' 

Pe 

Po 

q 

p,q,r,d,p 

P,4,r 


r ' 


r l 

S 


T l/2 


t 

V 

V 


W 


X 


0 


area -moment of the rudder about the hinge line 
mass of airplane, W/g, slugs 

period of damped natural frequency of airplane, sec 
rolling angular velocity factor, 2^, radians 
nozzle exit pressure, lb/sq ft, abs 
static pressure, lb/sq ft 
dynamic pressure, ipV^, lb/sq ft 

rate of change with time of cp, 0, ty, a, and p, 
respectively, radians/sec 

rate of change with time of p, q, and r, respectively, 
radians/sec^ 

yawing angular velocity factor, — , radians 
indicated yawing angular velocity, radians 
indicated yawing acceleration, radians/sec 2 
wing area, sq ft 

time required for absolute value of transient oscillation 
to damp to half amplitude, sec 

time, sec 

airspeed, ft /sec 

side velocity, ft/sec 

equivalent side velocity, (v{a ) , ft/sec 
weight of airplane, lb 

distance from center of gravity to transverse accelerometer 
(measured parallel to body X-axis), positive when forward 
of center of gravity, ft 

distance from center of gravity to sideslip vane (measured 
parallel to body X-axis), positive when forward of center 
of gravity, ft 
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distance from center of gravity to transverse accelerometer 
(measured perpendicular to tody X-axis) , positive when 
"below center of gravity, ft 

angle of attack of airplane, deg 

trim angle of attack for 1 g, deg 

corrected angle of sideslip,, deg or radians 

indicated angle of sideslip, radians 

total aileron deflection, positive when left aileron 
deflected down, deg 

elevator deflection, positive when elevator deflected 
down, deg 

rudder deflection, positive when rudder deflected to 
left, deg 

angle between referenced body X-axis and principal X-axis, 
positive when reference axis is above principal axis at 
the nose, deg 

ratio of actual damping to critical damping 

angle of pitch, positive when airplane nose pitches up, 
radians 

mass density of air, slugs /cu ft 
air -density ratio 

time conversion parameter, m/pVS, sec 
phase angle, deg 
damping angle, deg 

angle of roll, positive when right wing moves down, 
radians 

angle of yaw, positive when airplane turns to right, 
radians 
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The symbol | j | represents the absolute magnitude of a j quantity 
and is positive. When employed in an equation, the equation is considered 
to be a vector equation. 

The phase angle of a vector j relative to another vector k is 
indicated by the subscripts in The second subscript k is used as 

the reference . 


AIRPLANE 


The D- 558 -II research airplane (figs. 1 and 2) is a midwing airplane 
with a 55 ° swept wing and 40° swept tail surfaces as measured at 50 per- 
cent of the mean aerodynamic chord. It is equipped with plain flap-type 
unboosted control surfaces linked directly to the control wheel and 
rudder pedals. The D- 558 -II (145) airplane used for subsonic and tran- 
sonic investigations is powered by a J34-WE-40 turbojet engine and an 
LR8-RM-6 rocket engine. The D- 558 -II (l44) airplane used primarily for 
supersonic investigations is powered solely by an LR8-RM-6 rocket engine. 


In the D- 558 -II (145) airplane the turbojet engine exhausts out the 
bottom of the fuselage between the wing and the tail at an angle of 
approximately 8 ^ relative to the body X-axis. The airplane is powered 
by only the turbojet engine at Mach numbers below O. 85 . At Mach numbers 
above 0 . 85 , it is powered by both the turbojet and the rocket engines. 

The rocket engine in both airplanes has four nozzles which exhaust 
out the rear of the fuselage, with the thrust line of the combined nozzles 
coincident with the fuselage center line. These nozzles may be fired in 
any desired combinations. For some of the test flights the all-rocket 
airplane was equipped with nozzle extensions designed to expand the 
exhaust gases to ambient pressure at 28,000 feet. Figure 3 is a photo- 
graph of these nozzle extensions. 

The physical characteristics of the D-558-II airplane are presented 
in table I. 

The weight of the airplane at the time of the pulse maneuvers was 
between 10,000 and 12,000 pounds. The centers of gravity and moments of 
inertia for these weight values are listed in table II. 

The values of Iy were determined by adding the analytically deter- 
mined contributions of the fuel to the empty weight values of Iy deter - 

by ground— conducted oscillation techniques • The rather uncertain 
values of Iy determined by ground- conducted oscillations made it advis- 
able to use the manufacturer's estimate of I x for the design weight 
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condition, then analytically correct this value of I x for fuel con- 
sumed. The value of I z was determined by assuming I z = Ix + Iy 

INSTRUMENTATION AND INSTRUMENT ACCURACY 


Standard NACA instruments were used to record airspeed, altitude, 
angular velocities and accelerations, normal acceleration, transverse 
acceleration, angles of attack and sideslip, and rudder, aileron, eleva- 
tor, and stabilizer positions. The airspeed, altitude, and angles of 
attack and sideslip were sensed on the nose boom. All records were 
synchronized at 0.1-second intervals by a common timing circuit. 

The turn meters used to measure the angular velocities and accelera- 
tions were referenced to the body system of axes of the airplane and are 
considered accurate to within ±1.0 percent of scale range. The pitch 
turn meter had 2.1° uptilt relative to the X-axis; this error in instal- 
lation was not significant. The roll turn meter bad 0.5° uptilt relative 
to the X-axis; the yaw turn meter had 1.0° uptilt relative to the X-axis 
and 0.5° tilt to the left relative to the Y-axis. 

Indicated sideslip angles and angles of attack, measured by vane- 
type pickups, were corrected for roll and yaw rate, and pitch-rate effects, 
respectively. The pickups were magnetically damped and had dynamically 
flat frequency-response characteristics over the frequency range of the 
airplane. The pickups are statically accurate to ±0.2°. 

The ranges and dynamic characteristics for the pertinent instrumen- 
tation employed are: 


Function 

Range 

Undamped natural 
frequency, cps 

Damping 

ratio 

a, deg 

-20 to 40 

10.0 

0.70 

P, deg 

±50.0 

10.0 

.70 

p, radians /sec 

±3-5 

18.5 

.67 

q, radians/sec 

±1.0 

18.5 

.65 

r, radians /sec 

±1.0 

14.2 

.65 

a n , g units 

-1 to 8 

16.0 

.65 

a^, g units 

±1.0 

13.5 

.65 


All data employed in the analysis were corrected for instrument 
phase lag. Position corrections were applied by time-vector methods of 
analysis to indicated linear accelerometer readings. 
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Rudder, aileron, elevator, and stabilizer positions were measured by- 
standard control-position transmitters linked directly to the control sur- 
faces. The elevator position was measured referenced to the stabilizer. 

The transmitter -recorder system had a flat dynamic response over the fre- 
quency range of the control movements encountered and is considered 
accurate to within ±0.1°. 

The nose-boom installation for measuring the airspeed was calibrated 
by the MCA radar -phototheodolite method. The Mach numbers presented are 
considered accurate to ±0.01 at speeds below M = 0.90 and at speeds 
above M = 1.04; ±0.02 to ±0.03 in the Mach number range from 0.90 to 1.04. 


The test procedure consisted of recording the airplane response to 
abrupt elevator pulses in the longitudinal stability investigation and 
to abrupt rudder or aileron pulses in the lateral stability phase of the 
investigation. The maximum oscillatory sideslip was generally of the 
order of 3° or less for subsonic conditions and 2° or less for supersonic 
conditions. The difficulty of controlling the all-rocket airplane in 
supersonic flight (ref. l) made the investigation of the power effects, 
particularly on the lateral stability of this airplane, a lengthy process 
of repeated flights to obtain data which could be analyzed quantitatively. 
Initially, the pilot found it impossible to maintain neutral rudder after 
making supersonic power-on rudder -pulse maneuvers. It was necessary to 
devise a pin and slotted sector to enable him to hold the rudder in the 
neutral position following the pulse. This difficulty was not evident 
during the power-off supersonic flight. To eliminate any inadvertent 
aileron motions the control wheel was locked by a chain device for the 
supersonic maneuvers. 

The rocket engine chamber pressure was maintained at approximately 
245 psi gage. The nozzle exit pressure without nozzle extensions was 
l/l7 of the absolute chamber pressure; with nozzle extensions the nozzle 
exit pressure was 1/53 of- the absolute chamber pressure. The use of 
nozzle extensions provided a convenient means of changing the pressure 
ratio p e yT?Q* 


TESTS 


AMLYSIS OF FLIGHT DATA 


The longitudinal derivatives , 


a 



mined by using the following relations: 
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The lateral derivatives were obtained by employing the time-vector 
method outlined in reference 4. The results in the present paper, how- 
ever, are relative to the body axes, whereas reference 4 dealt with the 
stability axes. Figure 4 is a sample of the time-vector solution rela- 
tive to the body axes. The use of the body axes affected the vector 
diagram for the side-force equation (fig. 4(c)) by introducing an addi- 
tional vector |-2m 


There is considerable uncertainty regarding the inclination of the 
principal axis. Preliminary calculations indicated that an inclination 
of 1.5° down at the nose relative to the body X-axis provided reasonable 
correlation between flight and preliminary theoretical estimations of 
C n in the supersonic range. It is believed the inclination of the 

principal axis thus determined is known possibly to within 0.5°. However, 
this magnitude of uncertainty can result in poor accuracy of - C n ^ 

and still not appreciably affect the magnitude of C n ^, C^, and C^. 

Therefore, the lateral derivatives presented in this paper do not include 



The values of C 7 and C„ (fig. 5) required for the time-vector 

solution of the remaining derivatives were obtained from theoretical 
estimates . 


PRESENTATION OF RESULTS 


The results of the investigation are presented in two general parts: 
(l) the longitudinal stability characteristics, and (2) the lateral 
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stability characteristics. Inasmuch as particular emphasis is placed on 
power effects, the discussion, with the aid of figure 6, briefly reviews 
the mechanics of the power effects first. The pertinent results of the 
investigation are presented in the figures listed in the following 
tabulation: 

Figures 

Longitudinal stability characteristics: 

Time histories of elevator -induced oscillations 7 

Variation of trim and trim a with Mach number ... 8 

Period and damping characteristics 9 

Longitudinal stability derivatives 10 

Comparison of flight derivatives with wind-tunnel 

derivatives 

Lateral stability characteristics: 

Time histories comparing power-on and power-off 

lateral oscillations 12 

Period and damping characteristics 13 

Power effects on amplitude ratios and phase angles .... 14, 15 

Power influence on stability derivatives 16 

Comparison of lateral oscillatory characteristics 

with lateral damping criteria 17 

Time histories comparing power-on rudder -free and 

rudder -fixed oscillations 18 

Influence of free rudder on period and damping 

characteristics j power on 19 

Influence on stability derivatives of free rudder 

during powered flight 20 

Approximate variation of rudder hinge -moment parameter 

with Mach number at different pressure ratios 21 

In the appendix may be found: 

Vertical-tail areas used in estimating derivatives .... 22 

Comparison of flight derivatives with calculated and 
available wind-tunnel derivatives for power-off 

conditions 23 

DISCUSSION 

Before discussing the jet exhaust effects on the stability character- 
istics of the airplane, it is believed pertinent to consider some of the 
qualitative aspects of the shock-wave behavior at the juncture of the 
vertical tail and fuselage due to jet streams exhausting at the rear of 
the fuselage. On this basis it is desirable to review the principle of 
the shock -wave behavior. The following discussion is based primarily on 
reference 2. 
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When a body has a jet exhausting at a sufficiently high pressure 
ratio into a free stream in yawed supersonic flight, the stream flow 
cannot turn through the deflection angle formed by the body and the out- 
line of the jet exhaust. As a result, the boundary layer on the lee side 
of the body thickens and separates upstream of the jet exit (fig. 6). 

The external shock wave at the jet exit on the lee side of the body 
becomes a lambda shock with one leg forming forward in the lower pressure 
region and, in the absence of the vertical tail, the other leg remaining 
at the lip of the exit. On the upstream side of the body there is no 
separation of the boundary layer and a simple shock wave is present at 
the lip of the jet exit. 

The introduction of a vertical tail (or a horizontal tail) in this 
unsymmetrical flow field causes the lambda shock wave to move forward 
into the low-pressure region of the vertical tail; the shock wave on the 
other side of the vertical tail remains attached to the jet exit. Inas- 
much as the pressure behind the external shock wave is higher than in 
front of it, the forward movement of the lambda shock wave on the lee 
side of the vertical tail results in higher pressures on the vertical 
tail than would be experienced if the shock wave were not present. The 
result is a decrease in the stabilizing action of the vertical tail for 
rudder -locked conditions. If the configuration has some of the rudder 
surface within the field of the shock-wave action, the higher pressures 
behind the shock wave on the lee side of the rudder tend to turn the 
rudder into the free stream. With a free rudder, this shock-wave -induced 
turning of the rudder into the free stream provides a stabilizing action. 


Longitudinal Stability Characteristics 

Several time histories are presented in figure 7 and indicate no 
unusual longitudinal characteristics in the D-558-H airplane. The 
variation with Mach number of the level-flight lift coefficient 

and trim angle of attack oq, shown in figure 8 for the five test alti- 
tudes. indicated normal characteristics. The scatter of the C T test 

data in figure 8 is attributed to weight and altitude differences. The 
a,Q test points have been reduced to 1 g load -fact or conditions by applying 

the following correction: 


t£L = 


- !) 


The period and damping variation with Mach number and altitude 
(fig. 9) and the longitudinal derivatives (fig. 10) showed no power-on 
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jet exhaust effects for the nozzle -extension-on and nozzle -extension-off 
configurations which were used to obtain data for the pressure ratio 
(p e /Po) ran S e from about 2 to 15 . It appears likely that the horizontal- 

tail surface is located at a sufficient distance above the rocket nozzles 
and shock-wave formations to be out of the field of influence. Had the 
horizontal tail been located in the vicinity of the longitudinal refer- 
ence axis of the airplane with its trailing edge close to the edge of the 
body, it is conceivable that power effects would be evident. 

The period and damping ratio curves (fig. 9) show in general a 
normal variation with Mach number . Inasmuch as the derivatives shown in 
figure 11 do not show any altitude effects, the variations of the period 
and damping ratio with altitude are functions of air density for all 
practical purposes. Note, however, that in the Mach number range of O .85 
to 0.90 there are discontinuities which are reflected in the derivatives. 
The damping ratio curves (fig. 9) indicate a large decrease in damping 
ratio is experienced in the supersonic region. This decrease would be 
attributed primarily to the large negative increase in in the 

supersonic region over the values in the subsonic region (fig. 10 ) . 

The erratic behavior of the ^ curve (fig. 10) in the 

transonic Mach number range from 0.80 to about 0.92 is shown in refer- 
ence 5 to be a characteristic of wings of this type. 

A comparison of flight and wind-tunnel values of and is 

presented in figure 11. The flight -determined values of the derivatives 
show generally good agreement with the wind-tunnel data of references 6 
to 8 . 


Lateral Stability Characteristics 

Rudder fixed.- Time histories of representative rudder -pulse maneuvers 
for power-on and power-off conditions with the rudder fixed during the 
transient portion of the maneuvers are shown in figure 12. The destabi- 
lizing influence of power on the period is very evident. In any consider- 
ation of the power effects on the airplane, the cylinders fired and the 
pressure ratio P e /P 0 of P rime importance. Although the data 

obtained were not of the desired quantity, they did provide a good insight 
into the effect of cylinders fired and the influence of pressure ratios on 
the stability. It should be kept in mind that the following results of 
the present investigation are based on data wherein the maximum supersonic 
oscillatory sideslip angle was of the order of 2° or less. References 2 
and 5 point out that the magnitude of the jet exhaust effects (for rudder- 
free as well as rudder -fixed conditions) is influenced by the magnitude 
of the sideslip angle as well as pressure ratio and Mach number. 
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A study of the period curves in figure 13 for an altitude of approx- 
imately 60,000 feet shows that when number 1 cylinder was firing, there 
was no appreciable change in the period when the pressure ratio was 
increased from 4.7 to 15*2 at M = 1.4 and from 4.4 to 13 in the vicinity 
of M = I. 5 . When number 1 cylinder was eliminated from the firing com- 
binations as at a Mach number of 1.49 and I .58 (refer to period plot in 
fig. 13 ), the destabilizing influence of power decreased very markedly. 
Since the number 1 cylinder is in a dominant position for power effects 
(fig. 3 ), it would be the first cylinder to bring about a flow separation 
on the lee side of the vertical tail. Conversely, it would generally be 
supposed that a change in pressure ratio from about 4.5 bo about 14 would 
result in larger destabilizing influences than those shown. Had high 
pres sure -ratio data been available at M = 1.28 for comparison with low 
pressure-ratio data at this Mach number, greater changes in period might 
have been evidenced than shown at M = 1.4. This conjecture is based on 
limited data obtained at an altitude of 45,000 feet and will be discussed 
in a later section. 

A slightly favorable increase in ^ as a result- of power 

is indicated in the vicinity of a Mach number of 1.3 (fig. 13) J however, 
at the high test Mach number of I .56 the power effect appears to be 
detrimental. 


All the amplitude ratios and phase angles (figs. 14 and 15 ) show 
power effects which increase with increasing Mach number. The reduction 

in resulting from power (fig. 14) signifies a corresponding reduc- 

■f>l . .. |cp| 


(fig. l 6 ). The increase in 


tion in C-u \~-o- — /- — 

K I 

2.5 times the power-off value of 6.7 at M = 1.5 
in regard to handling qualities . 


which is approximately 
is usually undesirable 


The influence of jet exhaust effects on the lateral stability deriva- 
tives in which the number 1 cylinder was included in the firing combina- 
tion is shown in figure 16 . The seemingly moderate reduction in Cy^ 

at higher supersonic Mach numbers due to the power effects is evidenced 
as a serious reduction in C n ^ and a significant reduction in C^. For 

example, at M = 1.5 the decrease of 0.10 in Cy^ due to a loss in the 

vertical -tail effectiveness as a result of power effects, signifies a 
decrease in vertical -tail contribution to C n ^ from 0.34 to 0.27. Inas- 
much as the wing-fuselage combination provides a destabilizing contribu- 
tion to C n p of about - 0 . 19 , the seemingly small loss in vertical -tail 

contribution is very significant. 
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It will be noticed that the adverse effects of the jet exhaust on 
the stability derivatives tend to level off with increasing supersonic 
Mach number. It is quite possible that the power-on and power-off 
characteristics of the airplane may become very similar at some higher 
Mach number. 

No attempt is made at this time to consider the variation with Mach 
number of the lateral derivatives for power-off conditions. Consideration 
of the derivatives will be made later in sections dealing with the compar- 
ison of the derivatives with available wind-tunnel data and calculated 
derivatives . 

As a matter of interest the Dutch roll characteristics of the air- 
plane for power-off and power-on conditions of supersonic flight at an 
altitude of 60,000 feet are compared in figure 17 with appropriate phases 
of the damping criteria of reference 9 • The figure shows that neither 
power condition would satisfy boundary A at this altitude; however, if 
the airplane were considered to be equipped with artificial stabilization 
devices which were temporarily inoperative, the lower boundary shown in 
figure 17 for this condition would have been generally satisfied at an 
altitude of 60,000 feet and over the Mach number range investigated. The 
pilot considered the airplane easy to handle during power-off supersonic 
flight at an altitude of 60,000 feet and, even though the addition of 
power resulted in some deterioration of its handling qualities, it was 
still controllable. 

Rudder free.- As was explained in the "Tests" section, the pilot 
was unable to hold the rudder in a steady neutral position by using only 
the rudder pedals during the transient phase of a pulse maneuver in power- 
on supersonic flight. It was noted that this rudder -free, power-on 
lateral period of oscillation in the supersonic region was shorter than 
the rudder -locked, power-on lateral periods, thus suggesting an improve- 
ment in lateral stability during .oscillatory, power-on, supersonic flight 
by having the rudder free. This is shown in time histories presented in 
figure l8. 

The rudder -free, power-on data, obtained at 45,000 feet, are rather 
meager; however, the quality of the few data points utilized is considered 
to be good. The rudder -free, power-on stability characteristics based 
on these points are shown in figures 19 and 20. Included in these figures 
are several rudder -fixed points for both power conditions. Also included, 
for purposes of comparing trends and levels, are previously presented 
rudder -fixed curves for an altitude of 60,000 feet for both power 
conditions. 

Figure 19 and the derivative plots in figure 20 show that at a Mach 
number of 1.4l the rudder -free, power-on stability at the low pressure 
ratio of 2.2 with all cylinders firing was similar to the power-off 
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rudder -fixed stability. This signifies an improvement in rudder -free 
stability due to power effects, inasmuch as figure 21 (reproduced from 
ref. l) shows that the rudder tends to float with the relative wind during 
power-off flight at this Mach number. Thus the favorable effect of power 
on the rudder -free stability counteracts the adverse effect of power on 
rudder -fixed stability which was evident in figures 19 and 20 at a Mach 
number of 1.37 for the same power conditions. No rudder -free data were 
available for higher pressure ratios; however, on the basis of figure 21 
an increase in pressure ratio for rudder -free flight conditions would 
have resulted in further stabilizing effects up to a Mach number above 
1.6 (depending on pressure ratio), with the maximum influence occurring 
at a Mach number of about 1.4. 

It is evident that while the power effects cause a maximum positive 
floating tendency of the rudder at a Mach number of about 1.4 (fig. 21), 
the vertical tail as a whole (rudder -fixed conditions) is still experi- 
encing decreasing effectiveness. Although the inclination of the lambda 
shock wave tends to become normal to the surface at low supersonic Mach 
numbers, the shock wave is well to the rear and weak. As the Mach number 
increases to some higher value, it is believed that the observed power 
effects indicate that the increased strength and forward movement of the 
lambda shock wave more than compensate for the increased inclination of 
the shock wave. 

It was conjectured in the section on rudder -fixed stability that a 
change in pressure ratio from about 4.5 to 14 would probably result in 
larger destabilizing influences at a lower Mach number (such as 1.28) 
rather than at about 1.5. This conjecture was made on the basis of the 
very limited rudder -fixed test points at an altitude of 45,000 feet and 
Mach numbers of 1.27 and 1.37 in figures 20 and 21. These test points 
showed that an increase in pressure ratio P e /Po at M = 1.21 over the 

value at M = 1.37 more than offset the improvement in lateral stability 
which would have been obtained by not firing the number 1 cylinder at a 
Mach number of 1.28. 


CONCLUSIONS 


A flight investigation using pulse techniques has been made to 
determine the longitudinal and lateral stability characteristics of the 
D_ 558 -II airplane with particular reference to the jet exhaust effects 
of the rocket engine. The oscillatory sideslip motions were generally 
of the order of 3° or less for subsonic tests and 2° or less for super- 
sonic tests. The pertinent results of this investigation are summarized 
as follows: 
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1. Any cylinder firing combination tested that included the top 
cylinder resulted in a comparable loss in directional stability. These 
effects were most severe at the highest test Mach number of approximately 
1.6. With only the two middle cylinders firing (horizontal plane), the 
power effects were small. 

2. At the higher supersonic Mach numbers the large adverse power 
effects on directional stability were insensitive to pressure ratios 
varying from 4.4 to 15.2 in the Mach number range of 1.55 1° 1*56. 

3. Power effects cause the rudder to float into the relative wind 
during power-on yawed flight and, as a result, tend to offset the 
destabilizing influence of the jet exhausts evidenced during rudder- 
fixed yawed flight. 

4. The longitudinal stability of the airplane is not influenced by 
the jet exhausts primarily because the horizontal tail appears to be out 
of the field of action of the jet exhaust effects. 


High-Speed Flight Station, 

National Advisory Committee for Aeronautics, 
Edwards , Calif . , J une 18, 1957 • 


18 


NACA RM H57G09 


APPENDIX 

COMPARISON OF FLIGHT DERIVATIVES WITH WIND-TUNNEL 
DATA AND CALCULATED DERIVATIVES 


Preliminary comparisons of the power-off flight derivatives with the 
calculations of reference 10 showed rather large discrepancies at super- 
sonic speeds. It appeared that the overprediction of C n ^ resulted from 

an excessive estimate of the tail contribution to the overall stability. 
Consequently, it was decided to reestimate the tail contribution using 
a different effective tail area for supersonic conditions from that used 
for- subsonic conditions (fig. 22). The tail area for the supersonic 
calculations is in accord with the manufacturer's use of tail area for 
supersonic calculations for this airplane. The effective aspect ratios 
for the subsonic and supersonic tail areas were determined from refer- 
ence 11 to be 1.53 and I.V7, respectively. 


Tail-off estimates of C v , CL , and C 7 were based on wind-tunnel 

v p p 

data obtained from references 10 and 12. The wing contributions to the 
dynamic stability derivatives were estimated from the methods of refer- 
ences 13 to 19. Subsonic horizontal-tail-interference effects on the 
static and dynamic stability derivatives were estimated from references 20 
to 22. Vertical -tail contributions to the static and dynamic stability 
derivatives were calculated using effective aspect ratios, calculated 
lift-curve slopes (refs. 13 to 19), and the equations of reference 21. 


A comparison of the flight derivatives with available wind-tunnel 
derivatives (refs. 10, 12, and 23) and calculated derivatives is shown 
in figure 23. Although the subsonic wind-tunnel data are for M = 0.l6, 
the flight data extend to a sufficiently low Mach number of O.67 to show 
the degree of correlation. The low subsonic value of the flight deriva- 
tives CL and Ci shows good agreement with the wind-tunnel data. 

P P 

The low subsonic value of Cy^ from flight data is of lower magnitude 

than tunnel data; the reason for this discrepancy has not been determined. 
The subsonic value of CL from flight data at M = O.67 is at an angle 

H 

of attack of approximately 2.9°, indicating flight values of to be 


of lower magnitude than wind-tunnel data (-0.108 as against 
wind tunnel) . 


-0.123 for 


It is quite possible that the subsonic flight values of 
somewhat low in magnitude when compared to wind-tunnel data. 



are 


The reason 
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, and Iv, as 

. 1 91 

well as limitations in the accuracy of experimental determination of — 


for this may lie in uncertainties in the true value of 

racy of experimental determ: 

^ (fig. 5) have a significant influence on 
from flight data in the subsonic region; this 


The calculated values of 
the determination of C 


'P 


may be appreciated from a study of the sample vector solution for the 
rolling derivatives (fig. 4(e)). If Cj had been assumed equal to 

would have been -0.115 instead of -0.106. Errors in the 
values of I v or 


zero, 

“H 

Tl 

X wx 

I P I 

to bring about an incremental change of about ± 0.006 in the value of Cj 


of the order of ±5 percent would be sufficient 


At the higher supersonic Mach number of 1.6l, the power-off flight 
derivatives show good correlation with available wind-tunnel data. Inas- 
much as the tunnel data of reference 12 were based on a model equipped 
with the original (small) vertical tail, the correction to the Cy^ and 

C n p data at M = 1.6l to account for the enlarged tail was estimated 

on the basis of the experimental data of reference 24. 

The calculated derivatives of reference 10 and the calculated deriva 
tives of this paper are compared with flight derivatives in figure 25 . 

In the subsonic region both sets of calculated derivatives show good 


agreement with one another and with flight data for C n ^ and C \ 
agreement with flight data is moderately poor for Cy . For the 




P 


;p 


calculations of this paper show the greatest discrepancy with flight. 
This may be attributed in part to uncertainty of the position of the 
actual center of pressure of the vertical tail. 

In the supersonic region, the calculated derivatives of this paper 


are in closer agreement with flight data for Cy 


P' 


C n , and 
n p' 


'P 


than 


the extrapolation of reference 10, but in slightly poorer agreement 
for C 7 . 
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TABLE I 

PHYSICAL CHARACTERISTICS OF THE DOUGLAS D- 558 -II AIRPLANE 

Wing: 

Root airfoil section (normal to 50-percent chord of 

unswept panel) MCA 65 -010 

Tip airfoil section (normal to 50 -percent chord of 

unswept panel) MCA 63p-012 

Total area, sq ft 175*0 

Span, ft 25.O 

Mean aerodynamic chord, in 87 . 501 

Root chord (parallel to plane of symmetry), in 108.51 

Extended tip chord (parallel to plane of 

symmetry), in 6 l.l 8 

Taper ratio O .565 

Aspect ratio 5*570 

Sweep at 50-percent chord of unswept panel, deg .... 55*0 

Sweep of leading edge, deg 58*8 

Incidence of fuselage center line, deg 5*0 

Dihedral, deg -5*0 

Geometric twist, deg 0 

Total aileron area (rearward of hinge line), sq ft . . . 9*8 

Aileron travel (each), deg ±15 

Total flap area, sq ft 12 . 58 

Flap travel, deg 50 

Horizontal tail: 

Root airfoil section (normal to 50-percent chord of 

unswept panel) MCA 63 -010 

Tip airfoil section (normal to 30 -percent chord of 

unswept panel) MCA 63 -010 

Total area, sq ft 39*9 

Span, in 1^3*8 

Mean aerodynamic chord, in 41.75 

Root chord (parallel to plane of symmetry), in 53*6 

Extended tip chord (parallel to plane of 

symmetry), in 26.8 

Taper ratio 0.50 

Aspect ratio 3*59 

Sweep at 30-percent chord line of unswept panel, deg . . 40.0 

Dihedral, deg 0 

Elevator area, sq ft 9*^ 

Elevator travel, deg 

up 25 

Down 15 
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TABLE I.- Concluded. 

PHYSICAL CHARACTERISTICS OF THE DOUGLAS D-558-II AIRPLANE 


Stabilizer travel, deg 

Leading edge up 4 

Leading edge down 5 

Vertical tail: 

Airfoil section (normal to 30-percent chord of 

unswept panel) NACA 63 -010 

Effective area, (area above root chord), sq ft 3 6.6 

Height from fuselage reference line, in 98*0 

Root chord (parallel to fuselage center line), in. . . . 146.0 

Extended tip chord (parallel to fuselage center 

line ) , in 27-0 

Sweep angle at 30-percent chord of unswept 

panel, deg 49.0 

Rudder area (rearward of hinge line) , sq ft 6.15 

Rudder travel, deg ±25 

Fuselage: 

Length, ft 42.0 

Maximum diameter, in 60.0 

Fineness ratio 8.40 

Speed-retarder area, sq ft 5*25 

144 145 

Powerplant : 

Turbojet J34-WE-40 

Rocket LR8-RM-6 LR8-RM-6 

Airplane weight, lb: 

Full jet and rocket fuel 15,15± 

Full jet fuel 11,942 

Full rocket fuel 15,787 

No fuel 9,421 10,382 

Center -of -gravity locations, percent mean 
aerodynamic chord: 

Full jet and rocket fuel (gear up) 23-5 

Full jet fuel (gear up) 25.2 

Full rocket fuel (gear up) 24.6 

No fuel (gear up) 27*3 27. 0 

No fuel (gear down) 26.7 26.4 


4U 


NACA RM H5TG09 


25 


TABLE II 

MOMENTS OF INERTIA AND CENTERS OF GRAVITY FOR NOMINALLY 
EXTREME FLIGHT TEST WEIGHTS OF THE AIRPLANE 


Airplane 

D-558-n 

(1^4) 

D-558-II 

(145) 

We ight , lb 

10,000 

12,000 

10,000 

12,000 

T g"| 

2.920 

2.920 



J-y; ± 






55,500 

56,275 

51,500 

54,700 

I z , slug -ft 2 

56,220 

59,250 

— 

— 

Center of gravity, percent c . . . 

25 

25 

25 

25 
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Figure 


1.- Three -view drawing of the Douglas D-55 8 - H research airplane. 


504 



l-57-1650 

Figure 2 .- Photograph of the Douglas D-558-II research airplane without nozzle extensions. 
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E-1026 

Figure 3.- Photograph of the nozzle extensions on the Douglas D-558-II 

all-rocket research airplane. 
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Figure 4.- Typical sequence employed in the determination of lateral derivatives using flight 

data and the time vector. 
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Figure 4.- Concluded. 


NACA RM H57G09 


Cn 


P 


C l 


r 




Figure 5 .- Calculated values of and as functions of Mach number and altitude. 
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Figure 6.- Flow about a yawed jet exhausting at a high jet-to-stream 

pressure ratio into a free stream. Sketch duplicated from reference 2. 
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(a) Rocket and jet airplane; M = 0.6l; hp = 25, ^00 ft. 

Figure 7.- Time histories of longitudinal oscillations induced by an 

elevator pulse. 
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Figure 7 • - Concluded . 
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Figure 8.- Variation of 1 g lift coefficient and corresponding trim angle of attack with Mach 

number . 
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Figure 9*- Longitudinal period and damping characteristics of the D- 558 -II 
airplane as functions of Mach number and altitude. 
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Figure 10.- Variation of static and dynamic longitudinal stability 
derivatives of the D-558-H airplane with Mach number. 


v>l 

00 


Plight data 




Figure 11.- Comparison of longitudinal flight derivatives with available wind-tunnel derivatives 

for power on and power off. 
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(a) Power on, cylinders 1, 2, and k firing. 

Figure 12.- Time histories of lateral oscillations of the all-rocket 
airplane. Rudder held against slotted stop during transient 
oscillations . 
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(b) Power off. 


Figure 12.- Concluded. 
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Figure 15 •- Influence of power on the variation of lateral period and 
damping with Mach number. Control surfaces fixed. 
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Figure 14.- Influence of power on the variation of lateral amplitude 
ratios with Mach number. All-rocket airplane. Control surfaces 
fixed. 
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Figure 15 .- Influence of power on the variation of lateral phase angles 
with Mach number. All-rocket airplane. 
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Figure l6.- Influence of power on the variation of static and dynamic 
lateral stability derivatives with Mach number. Control surfaces 
fixed . 
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Figure IT.- Comparison of supersonic lateral oscillatory characteristics of the all-rocket air- ^ 

plane for power-on and power-off conditions with lateral damping criteria of reference 9* 
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(a) Rudder free; M = 1.41; h p = 44,500 ft. 

Figure 18.- Time histories of lateral oscillations of the all-rocket 
airplane for rudder -free and rudder -fixed conditions. All four 
rocket cylinders firing. 
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(b) Rudder fixed; M = 1-38; bp = 48,000 ft. 
Figure 18.- Concluded. 
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Figure 19.- Influence of free rudder during powered flight on the lateral 
period and damping of the all -rocket airplane. 


7U 


MCA RM H57G09 


1^9 



Rudder- fixed 






comparison curves Power 

h p , ft 

45,000 ft data 

Power 

Rudder 

• 

On, off 

30,000 

O 

Off 

Fixed 


Off> 


♦ 

On 

Fixed 


j 55 to 60,000 

4 - 

On 

Free 



M 

Figure 20.- Influence of free rudder during powered flight on the static 
and dynamic lateral stability derivatives. 
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Figure 21.- Approximate variation of rudder hinge-moment parameter with 
Mach number at different pressure ratios (reproduced from ref. 1). 
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Area = 56.6 sq ft 

(a) Subsonic calculations. 



(b) Supersonic calculations. 


Figure 22.- Vertical -tail-surface areas used in estimating values of 
calculated derivatives. Power off. 
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Figure 25 • - Comparison of lateral flight derivatives with calculated and 
available wind-tunnel derivatives for power-off conditions. Body axes. 


NACA - Langley Field, Va. 


